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^ ' Abstract 



Low energy supersymmetry with bilinear breaking of R-parity leads to a weak-scale 
seesaw mechanism for the atmospheric neutrino scale and a radiative mechanism for 
the solar neutrino scale. The model has striking implications for collider searches of 
supersymmetric particles. Assuming that the lightest SUSY particle is the lightest 
neutralino we demonstrate that (i) The neutralino decays inside the detector even 
for tiny neutrino masses, (ii) Measurements of the neutrino mixing angles lead to 
predictions for the ratios of various neutralino decay branching ratios implying an 
independent test of neutrino physics at future colliders, such as the Large Hadron 
Collider or a Linear Collider. We study the lightest neutralino decay branching ratio 
predictions taking into account present supersymmetric particle mass limits as well as 
restrictions coming from neutrino physics, with emphasis on the solar and atmospheric 
neutrino anomalies. 
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1 Introduction 



The simplest interpretation []1| of recent solar and atmospheric neutrino data [0, ^ ^ indicates 
that neutrinos are massive and, in contrast to the case of quarks, at least one and possibly 
two of the lepton mixing angles are large. The possibility of testing for these angles at high 
energy colliders seems very intriguing ||^, ^ especially in view of the new generation of 
colliders such as the LHC and a new high energy linear collider. 

One of the simplest ways to induce neutrino masses is if right-handed neutrinos, sin- 
glets under the SU{2) ® U{1) SM gauge group exist, as expected in a number of extended 
electroweak models and grand unified theories (GUTs) P, 0. In this case there are 
renormalizable neutrino Dirac masses similar to those of the charged fermions and, in addi- 
tion, potentially large Majorana mass terms for the right-handed neutrinos. This leads to a 
neutrino mass matrix of the form 



m 
m M 



(1) 



which has as eigenvalues 



M VM2 + 4rn2 
Ai,2 = y T ^ (2) 



where m = h'^ (H) with h"^ denoting a Dirac-type Yukawa coupling for neutrinos and (H) 
the vacuum expectation value (vev) of the SM Higgs field. It is also assumed that right 
handed neutrinos have a large mass term specified by M. This so-called seesaw-mechanism 
ll|| provides a general recipe to generate neutrino masses. For the simplest case m/M -C 1 



and just one generation one easily obtains that Ai ~ —im?/M while A2 — M -|- m? /M 
corresponds to a heavy right-handed state Q Typically the small neutrino masses required 
by the interpretation of present neutrino data correspond to values of M in a wide range 
above 10^ GeV or so. The parameters can be adjusted in such a way that they neatly 
explain neutrino experiments but they are far from being predictive. They mainly 'postdict' 
experimental data and lead to predictions, if any, which are confined to the domain of 
"neutrino" experiments, performed at underground installations, reactors or accelerators 
and possibly neutrino factories. They hardly imply any signatures that may be tested at 
high energy collider experiments such as the LHC. 

An alternative possibility for providing neutrino masses exists in which M is of order mz 
and that m is thus rather small. It is based on the idea of R-parity violation as origin of 
^For detailed results on diagonalization of seesaw neutrino mass matrices see WM- 



neutrino mass and mixing [1^, These models can have many imphcations for gauge 
and Yukawa unification, low-energy physics |T^, |T^, as well as a variety of implications 
for future collider experiments at high energies |1^, |18[. The simplest model in which 
this is realized is an extension of the Minimal Supersymmetric Standard Model (MSSM) |jl9[ 
with bilinear R-parity breaking terms in the superpotential [[14|, |15[] . 

These bilinear terms lead to non-vanishing vevs for the sneutrinos which in turn induce 
a mixing between neutrinos and neutralinos p, 0. This leads to an effective neutrino 



mass matrix which is projective [g, |2y] implying that only one neutrino receives mass at the 
tree-level. This provides for a way to account for the atmospheric neutrino problem. For 
the explanation of the solar neutrino puzzle one has to perform a 1-loop calculation of the 



neutrino/neutralino mass matrix in order to break the projectivity of the mass matrix P, pi 



The net effect is a hybrid scheme in which the atmospheric scale is generated by a weak-scale 
seesaw mechanism characterized by a mass scale is of order 10^ GeV, while the solar neutrino 
scale arises from genuine loop corrections ^ . 

In the bilinear Ijip model not only the neutrino masses but also the neutrino mixing angles 
are predicted in terms of the three underlying parameters |21| . These same parameters 
also determine the decay properties of the lightest supersymmetric particle (LSP) which we 
assume is the lightest neutralino. This implies the existence of simple correlations between 
neutrino mixing angles and neutralino decay branching ratios, as already partly observed 
1^. Note however that so far the literature [0 has focused mainly on qualitative statements 
based only on a tree- level approximation of the neutrino mass and of the neutrino/neutralino 
couplings. First of all this is not always reliable for extracting information on the atmospheric 
angle. In addition, it is totally unsuitable for making any determination of solar neutrino 
oscillation parameters from collider physics. 

In the present paper we present a quantitative approach to this problem which takes 
into account the complete 1-loop calculation of neutrino/neutralino masses and couplings. 
This is required to make reliable neutrino mass and mixing angle predictions. Moreover such 
complete treatment is necessary in order to establish reliable correlations between neutrino 
mixing angles and neutralino decay branching ratios. This also includes the solar mixing 
angle which has not yet been considered in this context in the literature. We discuss the most 
important correlations between neutrino physics and neutralino physics in detail. In addition 
we give an overview of all restrictions to the branching ratios of the lightes neutralino. 
Moreover, we have taken into account also LSP decays via real and virtual Higgs bosons, 
which have been so far neglected in the literature. Their contribution can be more important 
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than the one of the Z-boson if the hghtest neutrahno is mainly gaugino-hke as preferred by 
SUGRA scenarios and if third-generation fermions are present in the final state. This has 
some remarkable implications as we are going to demonstrate. 

Notice that we have ignored the results of the LSND experiment ||22| which would point 
to the existence of four neutrinos in Nature |23| one of which is sterile. Actually it is simple 
to extend our Ifip model in order to include a light SU{2) ®U{\) singlet superfield pi| . 
The fermion present in this superfield is the sterile neutrino, which combines with one linear 
combination of Vf. — v^ — to form a Dirac pair whose mass accounts for the LSND anomaly. 
On the other hand the sterile neutrino scalar partner can trigger the spontaneous violation of 
R-parity, thereby inducing the necessary mass splittings to fit also the solar and atmospheric 
neutrino data. This way the model can explain all neutrino oscillation data leading to 
four predictions for the neutrino oscillation parameters. However, for simplicity and for 
definiteness we will focus here on the simpler scenario with only the standard three light 
neutrinos. 

The paper is organized as follows: In Sec. we present the model and discuss approxima- 
tive formulas for some R-parity violation couplings. In Sec. |^ we discuss the phenomenology 
of the lightest neutralino at future accelerator experiments and in Sec. ^ we discuss in detail 
the relationship between neutrino-physics and neutralino-physics in this model. This in- 
cludes predictions for neutralino decays before SUSY is discovered and various cross-checks 
after the SUSY spectrum is known. In Sec. ^ we draw our conclusions. In Appendix ^ we 
give some formulas for approximate diagonalization of the generalized Higgs matrices. 



2 The model 

R-parity conservation is an ad hoc assumption in the MSSM and Ijip may arise either as 
unification remnant or through SIJ{7)®U{X) doublet left sneutrino vacuum expectation 



values (vevs) (i/j) |13|. Preferably we break R-parity spontaneously through singlet right 



sneutrino vevs, either by gauging L-number, in which case there is an additional Z |^ or 



within the SU{2) ® U{1) scheme, in which case there is a physical majoron In order to 
comply with LEP data on Z width we must assume that the violation of R-parity is driven 
by SU{2)®U{1) singlet sneutrino vevs p8[ since in this case the majoron has a suppressed 
coupling to the Z. Spontaneous R-parity violation may lead to a successful electro weak 
baryogenesis PSI . 



As long as the breaking of R-parity is spontaneous then only bilinear terms arise 
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in the effective theory below the Ijip violation scale. Bilinear R-parity violation may 
also be assumed ab initio as the fundamental theory. For example, it may be the only 
violation permitted by higher Abelian flavour symmetries . Moreover the bilinear model 
provides a theoretically self-consistent scheme in the sense that trilinear implies, by 
renormalization group effects, that also bilinear is present, but not conversely. 

The simplest Ifip model (we call it Ifip MSSM) is characterized by three independent 
parameters in addition to those specifying the minimal MSSM model. Using the conventions 



of refs. [|T9|, ^ the model is specified by the following superpotential [|T4 
W 



(3) 



where the couplings hjj, Hd and He are 3x3 Yukawa matrices and fi and are parameters 
with units of mass. The second bilinear term in eq. (^) includes R-Parity and lepton number 
violation in three generations. 

Similary the soft supersymmetry breaking terms are obtained by adding the correspond- 
ing R-parity breaking bilinear terms to the supersymmetry breaking Lagrangian of the 



MSSM. For the explicit form of these terms we refer to ||2T|. The important point for the 
following discussion is that beside the Higgs bosons also the sneutrinos aquire non-zero vevs, 
which we denote be vd =< >, vu =< >, Vi =< i>e >, ^2 =< > and ^3 =< >■ 
Note that the W boson acquires a mass = g^v'^/A, where = t>|, + + t>f + 1>| + 1>| ~ 
(246 GeV)2. Like in the MSSM we define tan/5 = vu/vd- 

2.1 Neutralino Mass matrix 

In the following we discuss the tree level structure of neutrino masses and mixings as 
needed for the following discussion. A complete discussion of the 1-loop mass matrix 

In the basis ip^'^ = 



and the other mass matrices in this model can be found in 21 



[—iX'i H^, H^, z/g, i/^, Ut-) the neutral fermion mass matrix M^v is given by 



M 



N 



m 



(4) 



where 





Ml 





-\g'vD 


\g'vu 







M2 


\gvD 


-\gvu 


- y'vD 













-\gvu 


-/^ 






(5) 
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is the standard MSSM neutralino mass matrix and 





-\g'vi 


\gvi 





ei 


m = 


- 


\gv2 



















(6) 



(7) 



characterizes the breaking of R-parity. The mass matrix Mtv is diagonahzed by 

A/'*MivA/'"^ = diag(m^o,m,,J 

where = 1, ■ ■ ■ , 4) for the neutrahnos, and (j = 1, ■ ■ ■ , 3) for the neutrinos. 

We are interested in the case where the neutrino mass which is chosen at tree level is small, 
since it will be determined in order to account for the atmospheric neutrino anomaly. Under 
this assumption one can perform a perturbative diagonalization of the neutrino/neutralino 
mass matrix |T2|, by defining ||20| 

e = m ■ M-J (8) 

If the elements of this matrix satisfy <^ 1 Vij then one can use it as expansion parameter 
in order to find an approximate solution for the mixing matrix M . Explicitly we have 

g'M2^l 



^4 



2det{M^o 
gM.fi 
2det{M^o 



■A,: 



-A,, 



{g^Mi+g'^M2)vu 



4:det{M^o] 



-- + 

{g^M^+g''M2)vD 
4:det{My^o) 



A,; 



(9) 



where 



Ai = fiVi + vnei (10) 

are the "alignment" parameters. From eq. @ and eq. (|l^) one can see that ^ = in the 
MSSM limit where = and Vi = 0. In leading order in ^ the mixing matrix A/" is given by, 



Af*-- 




1 - 

-e 1 - 



The sub-matrices N and V^, diagonahze M.x° ''^eff 



N*M-^oN^ = diag( 



VjnieffVy = diag(0, 0, m,y), 
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(12) 
(13) 



where 



m 



eff 



4:det{Myo] 



(14) 



and 



m,. 



Trim, 



eff) 



\AeAr A^A, Al ) 
M,g^ + M^g'^ ^ 



4det{M^o) ■ ^^^^ 
Due to the projective nature of rrieff, only one neutrino acquires mass p. As a result 
one can rotate away one of the three angles in the matrix V^^tree, leading to |^| 



tree 



/ COS 9^ 



/ 1 

cos 6*23 — sin 623 
\ sin 6123 cos ^23 / 

where the mixing angles can be expressed in terms of the alignment vector A as 

A, 



13 -sin 6*13 \ 
1 
V sinei3 



(16) 



cos 6*13 / 



tan^ 



13 



tan^' 



23 



(A^ + A2) 

^_A, 
A^ ■ 



1 5 
2 



(17) 



As discussed in detail in |^1]] the inclusion of 1-loop corrections to the mass matrix lifts 
the degeneracy between these states. Only after including these corrections one obtains a 
meaningful angle in the 1 — 2 sector. Both features are required to account for the solar 
neutrino data. 



2.2 Approximate Formulas for neutralino couplings 

The set of Feynman diagrams involved in neutralino decays is indicated in Figs. |l|, and 
1^. Most of the relevant couplings involved have been given in appendix B of ref. and 
the remaining ones will be included in appendix B of the present paper. Even though these 
are sufficient for our calculation of neutralino production and decay properties, it is very 
useful to have approximate formulas for the neutralino couplings, since this allows some 
qualitative understanding of the correlations we are going to discuss. To achieve this we 
make use of the expansions for the neutralino mass matrix and also a corresponding one 
for the charginos as given in [pO[]p| . For this purpose we will confine ourselves to the tree- 



level neutralino/neutrino mass matrix and we refer to Sec. [4.1| for a short discussion of the 
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Note that one has to reverse the sign of the ei in pO| to be consistent with our present notation. 




Figure 1: Feynman graphs for the decay Xi ~^ h ^t- 



necessary changes once the 1-loop corrections to the mass matrix are included. However, we 
have used exact numerical diagonalizations and loop effects in the calculation of all resulting 
physical quantities presented in Sees. 3 and 4. 

One class of decays which is important are those involving a VT-boson, either virtual or 
real. The x\-W^-li couplings are approximatively given by: 



Ocnw 
Ril 



gh%VD 

2Det. 



gvDN^2 + M^Nu 



+g 



/i 

+ g^pvu) Nu + (/x + Ma) gvuNu 
2/iDet+ 



A,; 



(19) 



gK 

V2 



2Detr 



Det4- 2Det 



vu fg^Mi+g'^M2 , g^ \ ,r : VD{g^Mi+g'^M2] 



+ 



Ni3 + 



14 



(20) 



2 V 2Deto /iDet+y 4Deto 

Here Det+ and Deto denote the determinant of the MSSM chargino and neutralino mass 
matrix, respectively. Nij are the elements of the mixing matrix which diagonalizes the 
MSSM neutralino mass matrix. 
For the coupling Z-Xi-i^i we find 



/^nnz 



y^nnz _ r, 

'g {gM,N,2 - g'M2Nn) ^ g {g^M, + g'^M2) vpN^ j 



4cos6'vi/Detc 



4 cos 6' Deto 



(21) 
(22) 




b) 



X? 



d) 



X?(Pi 



f) 



XiiPi) 



CO . pO 






Figure 2: Generic Feynman graphs for semi-leptonic neutralino decays. 




b) 



X? 



CO . pO 




Figure 3: Generic Feynman graphs for invisible neutrahno decays. 



As aheady mentioned, the tree-level the states z/i and i>2 are not well defined. Therefore one 
has to consider the complete 1-loop mass matrix as it will be done in the numerical part 
in sections 3 and 4. However, as one cannot detect single neutrino flavours, in experiments 
one observes the decay of Xi ^ X + Ui summing over all neutrinos z/j. Therefore, for the Z- 
mediated decay the interesting quantity is J2i=i 3 \^lx°u- ^^d, in contrast to the individual 



Xi 



Z z/j decay rates, this only gets small radiative corrections. 



For the coupling Xr^i-{Si — h^) we get 



nnh 



OnnI 
111 



E^o sin a C2 C4 Cg 



-ee (A^ + kl) + Ae (e/,A^ + e^A,) 



+ 



/iv/AfTAf|A| 
-S2 [Al + a;) + A, {s^A^ + s^Ar) 



(23) 



9 



0^3"/^ = E.o f sin a C2 C4 + ) - D^, C4 cg | A| (25) 



with 



(S2,S4,S6) , (26) 

.and (27) 

(^f^Mi + fif'^Mz) {[co^avD + sinat;,7) (s-'A^n - gNx2) + 2/i(sinaA^i3 - cosaA^M)] 

8Deto 

(28) 



The quantities and Cj are parts of the mixing matrix for the neutral scalars, which is 
discussed in Appendix ^ . 

For the couphngs cIlz — di — one finds 

ndnL -6e (A^ + A^) + A^ (e^A^ + e.A.) 

/iVAf+AflAI 

O^-Jf = (30) 
QduL ^ ^.^ e^A -6^A. 

^/i^Af+A^ 

0^:2^ = (32) 

O^^" = (33) 

= %|A| (34) 

with G^o = [g'^Mi + ^'^M2)t;t//(4Deto) and i/^o = {'ig'^Mi + ^'^M2)/i/(6y2Deto). For the 
couphngs dRi — di — Vj one finds that 0^"^^ = Ofl^ and O^"^^ = 0^"^" as above but with 
iy^o ^ ^'^Ma/i/ (3 V2Deto). 

One can obtain the couphngs between lu-k-i^j by replacing ho He and H^o {g'^Mi + 
g''^M2) jj, /{2y/2Deto) in the above equations. For the case of Im-k-Uj one finds the couplings 
by replacing he He and H^o — > g '^M2/i/(\/2Deto). 

For the couplings uu — Ui ^ Vj one finds 

OZi = 0^^,t = OTli^ = O"^,,^ = (35) 
0^3^ = -/^^J/x?|A| (36) 
or/' = ^x?|A| (37) 
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with J^o = (c/^Mi + c/'^M2)fD/(2Deto) and J^o = {-^g'^Mi + c/'^M2)/i/(6y2Deto). For 
the couphngs Ujn-Ui-Uj one finds that O^"^^ = O^"^^ and O™^^ = O'^^ as above but with 



Jr.g 



-V2^''M2/i/(3Deto). 



For the couphngs Uj — — k one finds 

^Lkl, 



vkk TDU 

■' ^ fi 2fi Det_|_ 



(3^ 



RL 



V2DetH 



+ 



+ h'h^M2R% 



+ 



M2 
1 + — 



A,, 



Det, 



(39) 



For the couphngs dj — Uk — h one finds 



^LkL 



^Rkli 



h%hfvDRji 

y2Det+ 



. . q^vn / , M2 A A,- 



fi J Det^ 



.kkr^d^^^ 9'R^i^9j!!>uSA 



Det- 



(40) 
(41) 



In eq. (^) - (|41|) we have assumed that there is no generation mixing between the squarks 
implying that j = 1,2. 

Data from reactor experiments indicate that the mixing element Ues must be small [Q. 
This implies that |Ae| <^ |A2,3|. In the limit Ae/A2,3 some of the above formulas simplify 
to 



Onnh 
111 



Onnh 
121 



r~\nnh 

•-^isi 



OdnL 
Liui 



OdnL 
Liu2 



OdnL 
Liu-i 



, C2C4 ce sm a eg 
-^x1 I + ^2 



E^o ( sin a C2 C4 cq- 



(^T^a - e„A^ sgA. - S4A 



/i|A, 



+ 



|Ar| 



E^o sin a C2 C4 Cg 



e^A^ + e^A^ S4A^ + SgA^ 



+ 



/iVAfTAf V^i + a: 



OdnR 
Riui 



/^dnR 
'-^Riui 



/^dnR 
^Rius 



Later on we wih also use the so-cahed sign-condition |2T[], defined by 

Au 

Cr A^ 



(42) 
(43) 

D^o C2 C4 CQsJkl + kl (44) 

(45) 
(46) 

(47) 



(4^ 
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Its origin can be traced back to the above eq. 



(^7\). Assuming ~ e,- as indicated by 



unification and |A^| ~ |At-| as required by the atmospheric neutrino problem one sees easily 
from the above equations that either the e part0 of the couplings to the second or the third 



neutrino state is very small depending on the relative sign between A^ and At-. If A, 



-A, 



loop 



v 



X K 



i',tree 



(49) 



one can show, after a lengthy calculation [Q, that the resulting effective neutrino mixing 
matrix is given by 

/ cos 9i2 — sin 9i2 \ 
sin 6*12 cos 6i2 
1/ 

with nearly unchanged ^13 and 623- In contrast, if this sign condition is not fulfilled the 613 
and 623 angles get large corrections. One sees therefore that if the sign condition is satisfied 
the atmospheric and solar neutrino features decouple: the atmospheric is mainly tree-level 
physics, while the solar neutrino anomaly is accounted for by genuine loop physics in a simple 
factorizable way. Thus the sign condition is helpful to get a better control on the parameters 
for the solar neutrino problem. 



3 Neutralino Production and Decays 

In this section we will discuss the production and the decay modes of the lightest neutralino 
Xi- In order to reduce the numbers of parameters we have performed the calculations in the 
framework of a minimal SUGRA version of bilinearly ^pSUSY model. Unless noted otherwise 
the parameters have been varied in the following ranges: M2 and from to 1 TeV, mo 
[0.2 TeV, 1.0 TeV], Ao/niQ and Bo/rriQ [-3,3] and tan/? [2.5,10], and for the parameters, 
[A^/y/AfTAfl = 0.4 - 2, e^/er = 0.8 - 1.25, |Ae/A^| = 0.025 - 2, e^/e^ = 0.015 - 2 and 
|A| = 0.05 — 0.2 GeV^. They were subsequently tested for consistency with the minimization 
(tadpole) conditions of the Higgs potential as well as for phenomenological constraints from 
supersjTiimetric particle searches. Moreover, they were checked to provide a solution to both 
solar and atmospheric neutrino problems. For the case of the solar neutrino anomaly we have 
accepted points which give either one of the large mixing angle solutions or the small mixing 
angle MSW solution. 

We have seen in eq. ( pTSD that the atmospheric scale is proportional |A|^/Det(m^o). As 

has been shown in [^ |2l| this statement remains valid after inclusion of 1-loop corrections 
■^The A parts lead only to a renormalization of the heaviest neutrino state whereas the e part gives mass 
to the lighter neutrinos. 
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Figure 4: a) Am^^^ and b) neutralino decay length. 
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Figure 5: Production cross section for the process cr(e+e XiXi) ^ function of m^o 
a Linear CoUider with 1 TeV c.m.s energy. ISR-corrections are included. 
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provided that |e|^/|A| < 1 implying that 1-loop corrections to the heaviest neutrino mass 
remain small. As we have seen in section ( p.2| ), most of the couplings are proportional to 
|A|/A/Det(m^o) and/or ej//i. Although |A|/(v^M2/i) has to be small in order to account for 
the atmospheric mass scale (see Fig. (^)) the previously discussed couplings are still large 
enough so that the neutralino decays inside the detector, as can be seen in Fig. (§b). 

In Fig. 1^ we show the cross section cr(e"^e^ XiXi) fb for ^/s = 1 TeV. Assuming now 
that an integrated luminosity of 1000 fb~^ per year can be achieved at a future linear collider 
(see and references therein) this implies that between 10^ to 10^ neutralino pairs can 
be directly produced per year. Due to the smallness of the R-parity violating couplings, 
most of the SUSY particles will decay according to the MSSM scheme implying that there 
will be many more neutralinos to study, namely from direct production as well as resulting 
from cascade decays of heavier SUSY particles. From this point of view the measurement 
of branching ratios as small as 10~^ should be feasible. As we will see in what follows this 
might be required in order to establish some of the correlations between neutrino mixing 
angles and the resulting neutralino decay observables, which is a characteristic feature of 
this class of models. 

In this model the neutralino can decay in the following ways 

ViVjVk (50) 

Viqq (51) 

^ ^^Ipk (52) 

^ Ifqq' (53) 

Pil (54) 

In the following we will discuss these possibilities in detail except Xi ^ ^il because its 
branching ratio is always below 10^^. 

In the following discussion we have always computed the complete three-body decay 
widths even in cases where m^o has been larger than one of the exchanged particle masses, 
so that two-body channels are open. This has turned out to be necessary because there are 
parameter combinations where the couplings to the lightest exchanged particle are 0(10) 
smaller than the coupling to one of the heavier particles, implying that the graph containing 
the heavy particle cannot be neglected with respect to the lighter particle contribution. An 
example is the case of Z-boson and S'j'-mediated gaugino-hke neutralino decays discussed 
later on. Here Si denotes the lightest neutral scalar. In addition we want to be sure not to 
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Figure 6: Invisible neutralino branching ratio summing over all neutrinos. 

miss possibly important interference effects as there are several graphs which contribute to 
a given process. A typical example is the process x? Ij where 26 contributions exist, 
as can be seen from the generic diagrams shown in Fig. |l|. 

The first important question to be answered is how large the invisible neutralino decay 
modes to neutrinos can be. This is important to ensure that sufficient many neutralino 
decays can be observed. As can be seen from Fig. ^ the invisible branching ratio never 
exceeds 10%. The main reason for this behaviour can be found in the fact that for the 
SUGRA motivated scenario under consideration the couplings of the lightest neutralino to 
the Z-boson are suppressed. This and the comparison with other couplings will be discussed 
in some detail later on. 

The mainly "visible" nature of the lightest neutralino decay, together with the short 
neutralino decay path discussed above, suggests the observability of neutralino-decay-induced 
events at collider experiments and this should stimulate dedicated detector studies. 

In Fig. we show the branching ratios for the decays into qqvi. Here we single out the 
6-quark (Fig. (0a)) and the c-quark (Fig. (0c)) because in these cases flavour detection is 
possible. One can clearly see that for m^o < 1.1 the decay into hbvi can be the dominant 
one. The reason is that the scalar contributions stemming from 5*^, and/or hk can be 
rather large. This can be understood with the help of eq. (p9|)-(^) where terms proportional 
to /i£)ej//i appear. This kind of terms is absent in the corresponding couplings for the u-type 
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a) Br(x2 ^ 66 ^i) b) Br(x? ^ T.q=u,d,s n) c) Br(x? ^ cc^i yi) 




m^o [GeV] m^o [GeV] m^o [GeV] 



Figure 7: Neutralino branching ratios for the decays into qqi'i final states summing over all 
neutrinos. 



a) Br(x? ^ e± E QQ') b) Br(x? ^ /x± E QQ') c) Br(x? ^ r± E ^^0 




m^o [GeV] m^o [GeV] m^^o [GeV] 

Figure 8: Neutralino branching ratios for the decays into l^q'q final states summing over all 
q'q combinations. 
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squarks implying that the branching ratio for ccvi is rather small as can be seen in Fig. (|^). 
One can see in Fig. (^) and Fig. (|^) a pronounced 'hole' around 80-100 GeV. It occurs 
because for m^o > my/ the W becomes on-shell implying a reduction for these decays. This 
is compensated as the Z becomes on-shell. 

The semi-leptonic branching ratios into charged leptons are shown in Fig. H. The decays 
into /i and r are particularly important because, as we will see in Sec. ^, they will give a 
measure of the atmospheric neutrino angle. Note that these branching ratios are larger than 
10~^ and in most cases larger than 10~^, implying that there should be sufficient statistics 
for investigations. In case of the e final state it might happen that one can only give an 



upper bound on this branching ratio. This is just a result of the reactor neutrino bound p3 



Note that due to the Majorana nature of the neutralino one expects in large regions of the 
parameter space several events with same sign di-leptons and four jets. 

In Fig. ^ the fully leptonic branching ratios are shown. One can clearly see a difference 
between the branching ratios into channels containing different charged leptons of the same 
flavour, i.e. t~t^ versus fi^fi^ and e~e~^. This difference is due to the importance of the 
Si state which corresponds mainly to the lightest Higgs boson of the MSSM. We have 
found that for gaugino-like Xi the couplings — Xi — are in general larger than the 
corresponding — Xi — couplings. This can be understood by inspecting the formulas 



given in eq. (|2l|) - (g^) in Sec. |2.2| , in particular the parts proportional to in eq. ( ^31) 
- (p5D. Other reasons for having "non- universal" t~t^ , A*"/^^ a-nd e~e'^ couplings are the 
graphs containing W or charged sleptons as exchanged particle (see Fig. p. From eq. (|19D 
one can see that the coupling 0|j^^ is proportional to h% implying that they only play a role 
if a r is present in the flnal state. 

Notice also that the largeness of the branching ratios for neutralino decays into lepton- 
flavour-violating channels can be simply understood from the importance of and 
contributions present in Fig. 11] 



4 Probing Neutrino mixing via Neutralino Decays 

In this section we will demonstrate that neutralino decay branching ratios are strongly 
correlated with neutrino mixing angles. We will consider two cases: 1) The situation before 
supersymmetry is discovered. In this case we demonstrate that neutrino physics implies 
^The charged scalars are a mixture of the charged Higgs bosons and the charged sleptons, and in particular 
the later are the important ones 
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Figure 9: Neutralino branching ratios for the decays into various lepton final states summing 
over all neutrinos. 
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^LZ\,aT)-pTox.l^LZ\ 



Figure 10: Approximated coupling O'j^^^pp^^^ using formula eq. (^) divided by the exact 
calculated coupling as a function of the exact calculated coupling. The bright (dark) points 
are for > (<)0. 



predictions for neutralino decays which will be tested at future colliders. 2) The situation 
when the spectrum is known to the 1% level or better as could, for example, be achieved at 
a future linear collider []35| , p7| . In this case our model allows for several consistency checks 
between neutrino physics (probed by underground and reactor experiments 0, ^) and 
neutralino physics. Moreover, some neutralino decay observables are sensitive to which of the 
possible solutions to the solar neutrino problem is the one realized, i.e. they can discriminate 
large angles solutions from the small angle MSW solution. 



4.1 Before SUSY is discovered 

Let us first consider the situation before SUSY is discovered. Before working out the predic- 
tions for neutralino decays we would like to point out a fact concerning the 1-loop corrected 



neutrino/neutralino mass matrix. It has been noticed in |21[ that the sign of the /i param- 
eter determines to some extent how large the absolute radiative corrections are (see Fig. 5 
of Q. The reason is that depending on this sign the interference between the 1-loop 
graphs containing gauginos and the 1-loop graphs containing Higgsinos is constructive or 
^The important information is the relative sign between /i and the gaugino mass parameters Mi^2- Since 
in pl| as well as here we assume that Mi^2 > then the absolute sign of n becomes relevant. 
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Figure 11: Testing the atmospheric angle. In case of the dark (bright) points fi < (>)0. In 
the second figure we have taken only those points with | sin 26*1! > 0.1. 



destructive. 

This fact has of course implications on whether the approximate couplings presented 
in Sec. ^]2| remain valid after the 1-loop corrections are taken into account 0. A typical 



example is shown in Fig. ^ where the approximated coupling Of^^^^^pp^^^ divided by the 
coupling O^gY as a function of Of^. One clearly sees that for /i < the tree-level result 
is a good approximation to within 20%, but for /i > it can be off by a factor up to 5 in 
some extreme cases where a constructive interference between gaugino and Higgsino loops 
takes place. We have checked that the same is true for the other couplings involving either 
the charged leptons and/or z/3. 



As can be seen from the discussion in Sec. p.2| the approximate formulas depend on the 
SUSY parameters, in particular on the parameters of the MSSM chargino/neutralino sector. 
However, one can see that the ratios of neutralino partial decay widths or of its branching 
ratios is rather insensitive to the MSSM parameters. As has been pointed out in |^ the 
atmospheric angle depends on the ratio of A^j/A^. This ratio (at tree level) can be obtained 
by taking the ratio Ofl^/Of^^. This leads immediately to the idea that the semi-leptonic 
branching ratios into fi^qq' and T^qq' should be related to the atmospheric angle. This is 
clearly demonstrated in Fig. O where we show the ratio of the corresponding branching 



®0f course the couplings involving vi and/or 1/2 are exceptional ones, as the angle between these states 
is only meaningful after performing 1-loop corrections are included. 
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Figure 12: Testing the atmospheric angle. In case of the dark (bright) points ^ < (>)0. In 
the second figure we have taken only those points with | sin2^^|;| > 0.1. 
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Figure 13: Testing the atmospheric angle. In case of the dark (bright) points ^ < (>)0. In 
the second figure we have taken only those points with | sin2^^| > 0.1. 
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Figure 14: Testing the Chooz angle. In case of the dark (bright) points fi < (>)0. In b) we 
have taken only those points with | sin2^^| > 0.1. 

ratios as a function of tan^ (^^atm)- One sees that present data imply that this ratio should 
be 0(1). In particular, the relative yield of muons and taus will specify whether or not the 
solution to the atmospheric neutrino anomaly occurs for parameter choices in the "normal" 
range or in the "dark-side", i.e. tan^( 



'atm. 



< 1 or tan2(^,i^) > 1 p 



The observed width of the band simply expresses the residual SUSY parameter depen- 
dence, which comes partly from the 1-loop calculated mass matrix and partly from the 
different contributions to these decays. If for some reason | sin 26'^ | > 0.1 the dependence on 
the parameters in the 1-loop calculation is considerably reduced because the sbottom/bottom 
loop dominates. This leads to a stronger correlation as seen in Fig. (pA]b). The fact that for 
fi > the band is wider is a consequence of the discussion in the previous paragraph. 

In Fig. |T2| and Fig. ^ we show two additional ratios which exhibit also a correlation with 
tan^^^i^: Br(x? e±/i^E.^.) / Mx'i ^ r±gg') and Br(x? / Br(x? ^ 

T^qq'). The (nearly) maximal mixing of atmospheric neutrinos implies that several other 
ratios of branching ratios are also fixed to within one order of magnitude, see Table 1 and 



Fig. |T7. 



In this model the so-called Chooz-angle is given by |Ae/AT-| |21] where we already have 
used the fact that the atmospheric data implies |A^| ^ |A^|. The same discussion as in the 
previous paragraph is valid. This leads automatically to the correlation between Bt{xi 
e^qq') /Bt{xi —>■ r^qq') and f/gg which is shown in Fig. 0. For f/^g < 0.01 the correlation 
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Figure 15: The solar mixing angle as a function of eg/e^. 

is less stringent because it implies that the tree level couplings have to be rather small and 
therefore loop corrections are more important. Note that existing reactor data give the 
constraint on U^^ < 0.05 at 90% CL [Q. This in turn implies an upper bound of ~ 0.2 on 
this ratio of branching ratios. 



The discussion of the solar angle is more involved. As illustrated in Fig. |T5| this angle is 
strongly correlated with ee/e^ ratio. In order to get information on the from neutralino 
decays one must take into account that, as already mentioned, the solar angle acquires a 
meaning only once the complete 1-loop corrections to the mass matrix have been included. 
For an easier understanding we focus on leptonic decays of the type Xi ~^ If^J^k with i ^ j 
which depend on the Xi-W-lj^i and W-kj-Uk couplings. The way a correlation appears is 
non-trivial. To understand it note that the couplings W-li-Uj depend on the neutrino mixing, 
since one must use the mass eigenstates for the calculation of the partial decay widths and not 
the electroweak eigenstates. In addition the enter via the Vj-S^-U and Xi-S^-k couplings. 
Remarkably, despite the non-trivial way the parameters enter here, one still has some 



residual correlation with ratios, as displayed in Fig. |16[ This figure shows that, although 
one does not get a strong correlation in this case, one can still derive lower and upper bounds 
depending on tan^(6'so/)- For the favored case [0] of the large mixing angle solution one finds 
that Bt{xi — ^ erz/j)/Br(x5 f^Tt^i) is determined to be one to within an order of magnitude. 



For the general bilinear model the spread in Fig. |T6| is due to the lack of knowledge of 



SUSY parameters. As will be shown in the next subsection, a much stronger correlation 
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Figure 16: Testing the solar angle. In case of the dark (bright) points /x < (>)0. In a) we 
have taken e^A^j/(e^A^) > 0. and in b) e^A^/(eT-AT-) < 

appears once the SUSY parameters get determined. 

In Table ^ we list upper and lower bounds on several ratios of branching ratios which are 
required by the consistency of the model. The values in the table are hardly dependent on 
the solution for the solar neutrino problem. 

The values in Table ^ can be viewed as important consistency checks of our model. 
However, one can also have observables which are able to discriminate between large and 



small angle solution of the solar neutrino problem. In Fig. we show how several ratios of 
neutralino decay branching ratios can be used to discriminate between large and small angle 
solution of the solar neutrino problem. The numbers in Fig. |1^ correspond to the following 
branching ratios: 1 . . . Br(ggi/j) / Br(e^r^z/j), 2 ...Br(66z/j) / Br(yU^r^z/j), 3 ...Br(66z/j) 
/ Br(r-r+z/i), 4 . . . Br(e±gg') / Br(/i±gg'). 5 . . . Br(e±gg') / Br(r±gg')> 6 ...Br(e±gg') / 
Br(eV^z/,), 7 . . . Br(/i±gg') / Br(eV^z/,), 8 . . . Br(/i±gg') / Bi{e^T^Ui), 9 . . . Br(r±gg') / 
Br(eV^z/i), 10 . . . Br(r±gg') / Br(e±r^z/i), 11 . . . Br(e±/i^z/i) / Br(e±r^z/i), 12 . . . Br(e±r^z/i) 
/ Br(/x"*"/i~z/j), and 13 . . . Br(/i^r^z/j) / Br(r+r^z/j). In Br(gg^^ z/j) we have summed over 
M, d, and s. Also for the case of z/,- we have summed over all neutrinos. 
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Table 1: Ratio of branching ratios as required by the consistency of the model. In 
Br(gg^. z/j) we have summed over u, d, and s. Also in case of Ui we have summed over all 
neutrinos. 



Ratio 


lower bound 


upper bound 
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Br(ggz/j) / Br(/i''=gg') 
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3.5 
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Br(ggz/j) / Br(e+e^z^j) 
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35 


Br(ggz/i) / Br^e^fi'^Ui) 
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Figure 17: Predicted ranges for the ratios of various branching ratios. The dark stripes are 
the ranges if one of the large mixing solutions (LMA, LOW or just-so) is realized in nature, 
the bright stripes are if SMA is realized in nature. The various ratios are given in the text. 
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Figure 18: Correlations between ia,T?Oatm and ratios of branching ratio for the parameter 
point specified in the text assuming that 10^ neutrahno decays have been measured. The 
bands correspond to an 1-cr error. 

4.2 After the SUSY spectrum is measured 

In the previous section we have discussed the predictions which can be estabhshed between 
neutrahno decay branching ratios and neutrino mixing angles before the first SUSY particle 
is discovered. Let us assume now that the entire spectrum has been measured with some 
precision, e.g. at a future Linear Collider As a typical example we discuss the point 

M2 = 120 GeV, /i = 500 GeV, tan /? = 5, setting all scalar mass parameters to 500 GeV, and 
also the A-parameter is assumed to be equal for all sfermions A = —500 GeV. Note, that 
we have taken /i positive to be conservative, as this corresponds to a 'worst-case' scenario." 
There are at least two parameters which need to be measured precisely: tan/5 and | sin2^^| 
because the 1-loop mass matrix is dominated by the sbottom/bottom loop if at least one of 
these parameters is large. 

In Fig. |l^ — Fig. ^ the same relationships as discussed above are displayed assuming 
that the particle spectrum and the corresponding mixing angles are known to the 1% level 
or better. In addition we have assumed that 10^ neutrahno decays have been identified and 
measured. Taking at the moment only the statistical error this translates to a relative error 
on the branching ratio Br(X) of the form 1/ A/lO^Br(X). It is clear from these figures that 
there exist excellent correlations between the ratio of various branchings and tan^ Oatm as well 
as the parameter f/^g probed in reactor experiments. For the solar angle we observe a strong 
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Figure 19: Correlation between f/^g and the ratio Br(e^gg')/Br(r''=gg') for the parameter 
point specified in the text assuming that 10^ neutrahno decays have been measured. The 
band corresponds to an 1-cr error. 
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Figure 20: Correlation between tan^ 6*30; and the ratio Br(e^r^ z/j)/Br(/i^r^ z/j) for 
the parameter point specified in the text assuming that 10^ neutralino decays have been 
measured. The band corresponds to an l-cr error. 
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dependence on tan^ Ogoi for the case of large mixing angle solutions (LMA, LOW or vacuum) 
of the solar neutrino problem. For the small mixing angle MSW solution, even though the 
dependence on tan^ Osoi becomes, unfortunately, rather weak, the ratio of branching ratios 
for Br(e^r^ v.i) /'Qi{^^t^ z/,) is predicted with good accuracy for any tan^ Osoi ^ 0.1. 



5 Conclusions 

Supersymmetry with broken R Parity provides a predictive hierarchical pattern for neutrino 
masses and mixings determined in terms of just three independent parameters, assuming 
CP to be conserved in the lepton sector. This can solve the solar and atmospheric neutrino 
anomalies in a way that allows leptonic mixing angles to be probed at high energy colliders, 
providing an independent determination of neutrino mixing. Taking into account data from 
atmospheric neutrino experiments we have derived specific predictions for neutralino decays. 



illustrated in Fig. O, Fig. 12 and Fig. 13. Probing the solar angle is more involved due 



to the intrinsic spread in SUSY parameters, see Fig. |I6[ However, we have demonstrated 
that, with about 10^ neutralino decays, and a determination of the spectrum of the theory 
to within 1% or better, there are very stringent correlations between solar neutrino-physics 



and neutralino-physics (see Fig. ^O]). We showed that several ratios of neutralino decay 
branching ratios can be used to discriminate between large and small angle solutions of the 
solar neutrino problem. Therefore, the hypothesis that bilinear R-parity violation is the 
origin of neutrino mass and mixing can be easily ruled out or confirmed at future collider 
experiments. This statement is actually more general, to the extent that the bilinear model 
is an effective theory of a model where R-parity is violated spontaneously. 
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A Approximate Diagonalization of scalar mass matrix 

Let us assume that the following simplifying conditions hold 

\€i€j\ <^ l-Bfcefel 

\Vi\ < Vd,Vu 

< Ia*I 

2 1 cos a {BiSi - gzVuVi) + sin a {gzVDVi - ^iSi) \ < 



(55) 

(56) 
(57) 



cos^ aAMrad + 9z {vu COS a — vd sin a)^ + Bqh (cot (3 cos a — tan /3 sin a)^ (vd/x — BiVu) 



(58) 



2| - sin q; (SjCj - gzVuVi) + cos a {gzVoVi - i^ei) \ < 



I sin^ aAMrad + gz {'Vu sin a + vd cos ck)^ + So/i (cot /3 sin a + tan /3 cos a)^ (w_d// — BiVu) \ 

(59) 

where gz = {g^+g'^) /4, AM^^^ = S^^^^/^ 157^2^2^ ^ g-^2 ^-^j^ g-^2 ^ ^ (v^+vj)) /{v^+ 
VD + vf + v^+v"^) and a is the mixing angle that diagonalizes the upper left 2x2 sub- 
matrix which corresponds to the usual Higgs mass matrix in the MSSM limit. Under these 
approximations the mixing matrix reads 

^ r'r.r'wr'^? sin rv rr,r^r,r rnfi rv .So .s^ .<?^ \ 



R 



C2C4^CQSma C2C4C6COSQ; S2 54 

C1C3C5COSQ; — CiCaCssino; Si S3 

— si cosa — S2 sina si sin a — S2 cos a C1C2 

—S3 cos a — S4 sin a S3 sin a — S4 cos a C3C4 

y — S5 cos q; — se sin a S5 sin a — se cos ck 



S5 



C5C6 y 



(60) 



with 



S2j = 



COS a {BiCi - gzvijv,) + sin a {gzVoVi - l^-^i) 



cos^ aAMrad + gz {vu COS q; — sin a)^ + i^oA* (cot /? cos a — tan /3 sin a)^ — ^ {vd/J' — BiVu) 



(61) 



S2i-1 = 



sin a {BiCi - gzVuVi) + cos a {gzVDVi - l^-^i) 



sin^ aAMrad + fl'z {vu sm.a-\-VB cos ct) + B^jj, (cot /3 sin a + tan /3 cos a) — ^ {vdI^ — B^vu) 



(62) 
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where i = 1, 2, 3 and q = 1 — sf/2. For the case where ^ (fD/W — BiVu) are equal for all i, 
the eiej part in the mixing between the sneutrinos becomes important. Therefore 



(63) 
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(64) 



On the other hand if ^ [vd^i ~ BiVu) = 
I J {vDfi - B2Vu) - ^ {vdI^ - Bsvu) I one has 



R' 



fl 

1 





y 
























B2Vu) Vo. i^DfJ' - B^vu) and 



\eiej\ < 







1 



/l 

1 





\^ -S7 














C8 








(65) 



-sg C7C8 y 



^ (i;i5/i - 5ii;t/) - ^ (t;fl/i - Esi^c/) 
^ (t;z5/i - 5iW[;) - ^ (t;B;U - 53WC/) 



(66) 
(67) 

(68) 

We have checked, that the eigenvalues obtained with these mixing matrices agree with 



those obtained in |29| in lowest order of the R-parity breaking parameters. 
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B Couplings 



Most of the couplings necessary for the calculation of the neutralino decays have already 
been given in where also the definition of the mixing matrices are given. The remaining 
couplings involve S^-Ui-di and qj-q[-lf- Neglecting the generation mixing among sfermions 
the corresponding Lagrangian is given by: 

+Sudi {as-,PL + bs-^Pn) + h.c. (69) 

with 

as-^ = hlRlt (70) 



h-^ = {KiTiRltr (71) 

= hlj'iRirV,; (72) 

cL^ = -9{R^irua + {h'h'r{R%ru.2 (73) 

CLk = h^D^Rim, (74) 

c%u. = -9{R%rva + {h'ij'nR%rv,2 (75) 
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